On a new definition of quantum entropy 
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It is proved here that, as a consequence of the unitary quantum evolution, the expectation value of 
a properly defined quantum entropy operator (as opposed to the non-evolving von Neumann entropy) 
can only increase during non adiabatic transformations and remains constant during adiabatic ones. 
Thus Clausius formulation of the second law is established as a theorem in quantum mechanics, 
in a way that is equivalent to the previously established formulation in terms of minimal work 
principle [A. E. Allahverdyan and T. M. Nieuwenhuizen, Phys. Rev. E 71, 046107 (2005)]. The 
corresponding Quantum Mechanical Principle of Entropy Increase is then illustrated with an exactly 
solvable example, namely the driven harmonic oscillator. Attention is paid to both microcanonical 
and canonical initial condition. The results are compared to their classical counterparts. 
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Allahverdyan and Niuwcnhuizen [1[ have recently 
pointed out that the formulation of the second law of 
thermodynamics in terms of Minimal Work Principle 
holds as a quantum mechanical theorem. That is, the 
unitary law of quantum mechanical evolution implies 
that the work W done by an external source on an iso- 
lated quantum system is larger than or equal to the value 
W corresponding to the quasi-static limit. This is sum- 
marized as: 



W > W 



(1) 



The result in {T]) is of great theoretical and practical rele- 
vance as it shows in a clear way how thermodynamic laws 
already exist at the microscopic quantum level. This is 
clearly indicated also by the works of Jarzynski on non- 
equilibrium free energy differences 0] ■ These results and 
related ones are now forming the basis of a new and ex- 
citing field of research with a great potential impact on 
nano-science, known as quantum thermodynamics [3[. 

In this Letter we will contribute to the field by show- 
ing theoretically and with examples that the Clausius 
Principle of Entropy Increase also holds as a theorem in 
quantum mechanics. This is accomplished by introduc- 
ing a new definition of quantum entropy operator, which 
is alternative to von Neumann's one 41: 



5(f) = \n(M(t) + 1/2) 



(2) 



where A/"(f ) is the time dependent quantum number op- 
erator: 



K 



#(*) = I>IM><MI 



(3) 



The reasons for this choice will be clarified later. 

For sake of clarity let us state Clusius Principle of En- 
tropy Increase here [f|: For every non- quasi- static pro- 
cess in a thermally isolated system which begins and ends 
in an equilibrium state, the entropy of the final state is 



greater than or equal to that of the initial state. For every 
quasi-static process in a thermally isolated system the en- 
tropy of the final state is equal to that of the initial state. 
Thus, we are going to prove and show with an exactly 
solvable example, that, for every non adiabatic trans- 
formation the expectation value of the quantum entropy 
operator of Eq. ^ cannot decrease. Further for every 
adiabatic transformation it remains constant. Note that 
von-Neumann entropy {—Trp In p) cannot satisfy such 
prescription because it always remains constant. Our 
result is slightly more general than Clausius principle as 
we relax the requirement that the final state be an equi- 
librium state. 

Consider an isolated non-degenerate quantum system 
which is initially in a state described by a density matrix 
p(U) = J2k=o Pk\k, ti)(k, ti\ which is diagonal over the 
Hamiltonian cigenstates basis {|fc,fj)}. Let the density 
matrix eigenvalues be ordered in a decreasing fashion: 



Prn > Pn if TU < U 



(4) 



At time ti a time dependent perturbation (not necessar- 
ily small) is switched on so that the Hamiltonian of the 
system can be expressed at any time t > f , as: 



K 



#(t) = $>*(*) IMXM 



(5) 



fc=0 



As time passes transitions will occur between the quan- 
tum states according to the transition probabilities: 

\a kn {t f )\ 2 = \(n,t f \U(U,t f )\k,t i )\ 2 (6) 

Where U(U,tf) is the unitary time evolution operator. 
As a consequence the density matrix evolves to some 
p(tf). Eq. (fTJ) has been proved under these assumption 

with work defined as W = Tr p(t f )H(t f ) - p(t t )H(ti) 

Here we are interested in the change in the expecta- 
tion value of the entropy operator: 



Sf-S^Tr p(t f )S(t f ) - p(ti)S(ti) 



(7) 



2 



Let us denote the probability that the system is found 
in state n at time tf, by p' nl where evidently p' n = 



T,k=oPk\akn{t f )\ 2 , then 



A" 



Sf-S i = J2(Pn-Pn)Mn+l/2) 

n=0 

Using the "summation by parts" rule Eq. 
comes 



K 



S f - ^ = In ~~~T y^iPn - Pn) 

m— ^ n— U 



(8) 



be- 



(9) 



Using the property of the coefficients la^l 2 of form- 
ing a doubly stochastic matrix (i.e., J2k=o \ a kn{tf)\ 2 = 
S„=o \ a kn(tf)\ 2 = 1 and \a kn (tf)\ 2 > 0), and the or- 
dering of probabilities j4|, it can be shown [l|, |4J that 



Hn=o(Pn - Pn) > 0, which, by noting that ln^±£ > 
in Eq. ([9]), implies 



5/ > 5 t . 



(10) 



This holds as a consequence of the laws of quantum me- 
chanics for every time dependent perturbation acted on a 
non-degenerate quantum mechanical system. As such it 
is an exact non- equilibrium result. In case the perturba- 
tion is adiabatic, the quantum adiabatic theorem would 
ensure that no transition will occur between states with 
different quantum number Q so that p' n = p n and con- 
sequently Sf = Si. Eq. (|10[) is the Quantum Mechanical 
Principle of Entropy Increase. Note that, like Clausius 
Principle, Eq. (JTDJ) in no ways implies that the expec- 
tation value of the entropy is monotonic increasing. Eq. 
(|10[) says only that after time ti the expectation value of 
the entropy operator will never be less than the initial 
value. This does not rule out the possibility that for two 
times ti > ia > U one might have Si < Sa- 

The reason for defining the entropy operator as in 
Eq. ([2]) can be understood by considering the classical mi- 
crocanonical entropy 



S cl (E,t) = In 



§(£,t) = / 

J H(q,p;t; 



<E 



d 3N qd 3N p 



(11) 

(12) 



that is the logarithm of the volume $ of classical phase 
space enclosed by the hyper-surface of energy E. Using 
the semi-classical viewpoint [H, pp. 23-24] the integral in 
Eq. (|12p counts the number of quantum states not above 
a certain energy e n = E. Since the levels are assumed to 
be non degenerate, this number is n + i, where we set by 
convention that the vacuum state counts as a half state. 
Thus one can construct the quantum version of Eq. (fT2|) 
and obtain Eq. By repeating all the calculations re- 
ported above in the classical case one can prove that also 



the classical expectation value of the classical entropy 
cannot decrease Q- This requires replacing the discrete 
quantity n + 1/2 with the enclosed volume and the 
probability p n with the probability density function P($) 
of having one member of the ensemble being on an orbit 
that encloses a volume $. Accordingly all sums over n 
will be replaced with integrals over d<£>. The coefficients 
|afc„| 2 have to be replaced by the "classical transition 
probability" A($, 8) that a representative point which 
is on an orbit that encloses a volume $ at time ti will be 
found on an orbit that at time tf encloses a volume 0. 
Thus, if the initial distribution is -P(3>), the final one is 
P'(9) = / d$A(6,$)P($) [1]. Liouville's Theorem im- 
plies that the "classical transition probability" A($, 0) is 
doubly stochastic, i.e., / A(G, = f A(Q,$)dQ = 1 

If one translates the quantum mechanical requirement 
of non-degeneracy (only one state per energy eigenvalue) 
with the classical requirement of ergodicity (only one tra- 
jectory per energy level) [3] , and the requirement in Eq. 
dU with the requirement that P($) be monotonic de- 
creasing, then one obtains Q (compare with Eq. ©): 



Sa 



d$(P'($) - P($)ln$ > (13) 



Let us illustrate these results with a practical example 
of an exactly solvable time-dependent quantum system. 
Consider the driven harmonic oscillator: 



H(x,p,t) = ^ + ywV + f{t)x 



(14) 



which can be solved analytically and is relevant in many 
applications such as polyatomic molecules in varying ex- 
ternal fields, colliding polyatomic molecules and elec- 
trons in crystals [ll]]. F° r simplicity let us set m = 1, 
u> = 1, ti = 0, tf = T, and consider the case when 
the perturbation is cyclic and switched on only during 
the time interval [0, T][ll|. Let us pick a representa- 
tive point from our initial ensemble, located at Xo,po 
at time U, and draw in phase space the circular orbit 
H(x,p, 0) = H(xo,po,0). Let ro be the radius of such 

r 2 

orbit, then will be the energy Eq of the system and 
the enclosed volume $ will be given by ^-E = Let 
us then choose the units in such a way that h = 1 . Then, 
in this specific problem, energy and enclosed volume co- 
incide. Using Laplace transforms one can easily solve 
the classical problem and obtain the energy E at time T, 
as: E = E + W + 2 v / iWcos(T - <^ - 0) or, using the 
enclosed volume notation (i.e., Eq = <I>, E = 0): 



= $ + W(T) + 2 ^§W{T) cos(T - cp - 6(T)) (15) 
where 

2W(T) - |/3(T)| 2 , 6(T) = arg/3(T) (16) 

f3(T) = f dtf(t)e iut (17) 
Jo 



3 



Thus, depending on the initial phase </j, a particle of 
initial enclosed volume $ will have different final en- 
closed volume 8 at time T. The probability of end- 
ing up enclosing a volume Q from a point that ini- 
tially encloses a volume $, then is easily calculated as 



A(6, <&) = !/ 



./<-) 

dip 



Using Eq. (TT5)) we get: 



A(6,$) = - 



7T (6-$- VF) 2 



(18) 



One can check that A(<d, $) is doubly stochastic. Us- 
ing Eq. (TT8|) one finds the average microcanonical final 
enclosed volume and its square deviation: 



(e>< 



<e% - (e) 2 $ = 2*W 



(19) 
(20) 



where the symbol (-) $ denotes average over the initial 
microcanonical distribution of energy Eq = $. Note, 
from Eq. (|19p . that W is actually the work done on the 
initial ensemble. The microcanonical expectation of final 
entropy is most conveniently calculated directly from Eq. 
(dU as [H: 

(ln9) $ = / \n{$ + W + cos <p)d(p 

Jo 

= ln[max($, W)] (21) 

Note that trivially one has (lnO) $ > ln$. That is, Clau- 
sius formulation is satisfied in this case even for an initial 
microcanonical distribution. According to Eq. p^|) . W 
is definite positive, so also the Minimal Work Principle 
(TT]), which for a cyclic process reads W > (Thomson 
formulation) is satisfied. Let us calculate the entropy 
change resulting from an initial canonical distribution: 



From Eq. (TIB"]) we get: 



(22) 



Sf-Sf = d$e-^*[(lne) $ -In*] 
Jo 



-0W I dxe- 0Wx lnx 



(23) 



which is evidently positive as expected. 

Let us now move to the quantum mechanical case. The 
transition probabilities | pressed in terms of 

Charlier polynomials C(m,n\W) (TH Il3j: 



i(m,n) 



C(m,n\W)= J2 



{-l) l m\n\ 



l\{m - l)\{n - l)\W l 



(24) 



\a nm (T)\ 2 = - 1 j [C(m,n|lU(T))] 2 (25) 

— if n] 



The average microcanonical quantum number and its 
square deviation are analogous to the classical formulas 
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FIG. 1: Microcanonical expectation of entropy as function of 
the initial quantum number n for W = 10. Solid line: classical 
case, Eq. (|2ip . Dots: quantum case, Eq. (|28p . the sum 
has been truncated at m = 1000. The quantum phenomena 
smooth out the sharp angle in the classical graph. 



PJ] (compare with Eqs. i|19|20p ): 



(m 2 ) n - (m); 



n + W 
2(n+ 1/2)W 



(26) 
(27) 



where the symbol (•) denotes average over the initial 
quantum mechanical microcanonical distribution pi = 
|n, 0)(n, 0|. The average microcanonical quantum expec- 
tation of entropy reads 



(ln(m + l/2))„ 



E 



a„ m (i)| 2 ln(m + l/2) 



(28) 



Unlike the average quantum number and square devia- 
tion this does not correspond exactly to the classical ex- 
pression. Figure [1] shows the quantum and classical aver- 
age microcanonical expectation of entropy. Note that the 
quantum effects smooth out the sharp angle that appears 
in the classical case. Note also that Clausius principle is 
satisfied for the microcanonical initial condition also in 
the quantum case, that is (ln(m + l/2)) n > ln(n + 1/2). 
Since any distribution is a superposition of microcanon- 
ical distributions, then Clausius formulation is satisfied 
here for any initial distribution, without restriction to 
decreasing ones This is true for the Minimal Work 
Principle as well (|26|) . Let us calculate the change in 
quantum expectation of entropy resulting from an initial 
canonical distribution: 



ft = (l-e-^e-""| n ,0)(n,0| 



(29) 



we have: 

Sf-Si = (1 -e- p )Y,cr Pn x (30) 

n 

[(ln(m + l/2)) n - ln(n + 1/2)] > 
Let us now study the actual change in expectation value 
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FIG. 2: Change in microcanonical expectation of entropy for 
the force in Eq. (|31[1 with L — 6 , n — 2, as a function of 
switching time T. Solid line: classical case, Eq. (|2ip . with 
$ = n + 1/2 = 5/2. Dashed line: quantum case, Eq. {28$, 
the summation is truncated at m = 1000. 




The Quantum Mechanical Principle of Entropy In- 
crease has been established here by replacing von- 
Neumann entropy with a new quantum entropy operator 
which is the natural quantum version of the classical mi- 
crocanonical entropy. The unitary character of the quan- 
tum time evolution is responsible for keeping the former 
always constant and making the latter grow. The va- 
lidity of the Quantum Mechanical Principle of Entropy 
Increase has been checked on the driven harmonic oscil- 
lator. Interestingly, we have noticed that the principle is 
satisfied in this example even in the microcanonical en- 
semble, thus showing that the requirement of decreasing 
probabilities is sufficient but not necessary. Further in- 
vestigations will aim at establishing under what general 
conditions the principle is satisfied in the microcanoni- 
cal case. Fundamental implications, equivalence with the 
Minimal Work Principle, and the mirror-image Principle 
of Entropy Decrease (valid in the negative temperature 
scenario), have been discussed elsewhere [H, E3]. These 
results contribute to our understanding of quantum ther- 
modynamics. The Quantum Mechanical Principle of En- 
tropy Increase could be tested experimentally, for exam- 
ple, with single molecule pulling experiments. 

An award from the Texas Section of the American 
Physical Society for the presentation of this work at their 
Fall '07 meeting is gratefully acknowledged. 



FIG. 3: Change in canonical expectation of entropy for the 
force in Eq. (|31[) with L = 6, f) = 2, as a function of switching 
time T. The canonical distribution is truncated at n = 100. 
Solid line: classical case, Eq. (|21|l . with $ = n + 1/2 = 
5/2. Dashed line: quantum case, Eq. (|28p . the summation is 
truncated at m = 1000. 



of the entropy for a specific shape of f(t). Let 



/(*) 



W(T) 



Lsin^j t€[0,T]; 

L 2 7r 2 r 2 (l+cosT) 
(tt 2 - T 2 ) 2 



otherwise (31) 
(32) 



Figure [2] shows the quantum and classical microcanon- 
ical entropy change resulting from the force (|3ip . as a 
function of switching time T. Note that as expected 
the change in entropy goes to zero in the adiabatic 
limit T — > oo. Note also a curious fact. The func- 
tion W(T) is bounded from above W{T) < cL 2 . Thus, 
for a given initial energy $ and amplitude L, if L is 
not sufficiently large, the change in microcanonical en- 
tropy remains zero for all values of switching time T\ 
This is evident in the classical case for which one has 
W(T) < $ =S> (ln6) $ = ln$,VT . 



L < 

Figure [3] shows the quantum and classical canonical 
entropy change. Note that, in all cases, in accordance 
with Eq. (|10p. the entropy change is non-negative. 
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